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Abstract 

We introduce a new approach to isolate the real roots of a square-free polynomial F = LjLo^< x ' w ' tn rea l 
coefficients. It is assumed that each coefficient of F can be approximated to any specified error bound. 
The presented method is exact, complete and deterministic. Due to its similarities to the Descartes method, 
we also consider it practical and easy to implement. Compared to previous approaches, our new method 
achieves a significantly better bit complexity. It is further shown that the hardness of isolating the real roots 
of F is exclusively determined by the geometry of the roots and not by the complexity or the size of the 
coefficients. For the special case where F has integer coefficients of maximal bitsize T, our bound on the 
bit complexity writes as <5(n 3 T 2 ) which improves the best bounds known for existing practical algorithms 
by a factor of n = degF. 

The crucial idea underlying the new approach is to run an approximate version of the Descartes method, 
where, in each subdivision step, we only consider approximations of the intermediate results to a certain 
precision. We give an upper bound on the maximal precision that is needed for isolating the roots of F. For 
integer polynomials, this bound is by a factor n lower than that of the precision needed when using exact 
arithmetic explaining the improved bound on the bit complexity. 

Key words: Root isolation, complexity bounds, bitstream coefficients, approximate coefficients 



1. Introduction 

Finding the roots of a univariate polynomial FfMJx] can be considered as the fundamental 
problem of computational algebra, and there exist numerous approaches dedicated to approx- 
imate the real roots of F. We mainly distinguish between purely numerical methods such as 
Newton iteration and exact and complete methods such as those based on Descartes' Rule of 
Signs or Sturm Sequences. The latter approaches apply to polynomials with rational coefficients 
and guarantee to compute a set of disjoint isolating intervals. That is, each of these intervals 
contains exactly one root and the union of all intervals covers all real roots of F. In this paper, we 
propose an algorithm which extends the Descartes method to arbitrary square-free polynomials 
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with real coefficient. Throughout the paper, 

n 

F(x) := J^Aj-x' £ R[x] (1.1) 

denotes a square-free polynomial of degree n > 2 with real coefficients A,, where |A„ | > 1 . We 
define X to be the minimal positive integer with max^o,. Jj4 < 2 T . It is assumed that each 
coefficient A, can be approximated to any specified precision and we refer to such coefficients as 
bitstream coefficients. The roots of F are denoted by §i, . . . ,§„ £C and D? := log (max, de- 
notes the corresponding logarithmic root bound. The separation a, := a(§,-,F) of is defined as 
the minimal distance of to any root £,j ^ the separation Cp of F is defined as the minimum 
of all ct(&,/), and E f := - £? =1 log a t . 

1.1. Main results and related work 

We present an exact and deterministic algorithm which computes isolating intervals I\,...,I m 
for the real roots of F. We further provide a detailed complexity analysis showing that our algo- 
rithm needs no more than 

0{n{Y. F + nT F f) = 0(n(L F +nx) z ) (1.2) 

bit operation^ and demands for approximations of the coefficients of F to <5(£f + riFp) bits 
after the binary point. Our results show that the complexity of isolating the real roots does not de- 
pend on whether the given polynomial has irrational, rational or integer coefficients. In fact, the 
hardness of isolating the roots of F is exclusively determined by the degree of F and the quan- 
tities Ff and Lp which only depend on the location of the roots of F. For a polynomial F with 
integer coefficients, the bound in (1.2) writes as (5(n 3 T 2 ) which improves the best bounds known 
for other practical methods such as the Descartes method [2,8,14,25,32], Sturm's method [10,23] 
or the continued fraction method [1,36,38,39] by a factor of n. To the best of our knowledge, this 
is the first time where it is shown that approximation leads to a better worst case complexity for 
real root isolation, a fact which has already been observed in experiments [17,32]. We consider 
this new result as an important step to further reduce the gap (with respect to worst case bit 
complexity) between practical and efficient algorithms for real root isolation and asymptotically 
fast methods for isolating all complex roots as proposed by Schonhage [35] and Pan [30,31] in 
the eighties and nineties. The latter methods achieve almost optimal complexity bounds <3(n 3 f) 
for the benchmark problem of isolating all complex roots but both methods lack evidence of be- 
ing efficient in practice; see [16] for an implementation of the splitting circle method within the 
Computer Algebra system Pari/GP. Due to its similarities to the Descartes method, we consider 
the proposed algorithm practical and easy to implement. The latter claim has already been proven 
by means of a recent implementation from A. Strzebonski and E. Tsigaridas [37] "in C as part of 
the core library of MATHEMATICAL 

The crucial idea underlying the presented method is to use an "approximate version" of the 
Descartes method. More precisely, we first consider a scaled polynomial f(x) := F(2 r+l x)/A„, 
where F is an integer approximation of Fp with Fp < F < Fp +41og«; see Section 2.2 and Ap- 
pendix 6.1. Then, all roots of / are contained within the disc of radius 1 /2 centered at the origin. 



O indicates that we omit polylogarithmic factors 



2 



Fig. 1.1. The above figure shows the recursion tree induced by the Descartes method when applied to the 
polynomial f(x) := 16V2x 2 — 8x + j (with roots z\ = 0.06... and Z2 = 0.29...). For each interval / = (a, b) 
in the subdivision process, we have to compute fi(x) = f(a+(b — ajx). For instance, for / = we have 

fl(x) = /(j + | ) = \/2x 2 + (2\/2 — 2)x+ \/2+ 5. The Bitstream solvers presented in [24,33] initially start with an 
approximation g of / to a certain number of bits; e.g., g(x) = '^jp ' x 2 — %x+ §77 approximates / to 10 bits after the 
binary point. Then, the Descartes method is applied to g, that is, for each interval / = (a,b), gi(x) = g(a + (b — a)x) is 
computed; e.g., g(\ + ^) = iff^r* 2 + tuft* — gffl- Gi ven that g is a sufficiently good approximation of/, it is shown 
that the roots of / can be isolated in this way. Our new approach follows a similar strategy, that is, we start with an 
approximation fi of // to a certain number p/ = p of bits. Then, we recursively compute approximations /; of /; 
to p; bits, where p; is updated in each step. In contrast to the previous method, the polynomials /; do not necessarily 
correspond to a specific initial approximation g of /. We illustrate this by means of the above example: We start with 
4W = TT? W + Trf which approximates f, (x) = f{-\+x) = 16V2x 2 - I6y/2x + 4 - 8.t + 4 + f +4^2 
to p; = 10 bits. Then, /; ( j ) and fi (j + j) are evaluated and the result is rounded to 9 bits after the binary point. The 
resulting polynomials are then approximations of (x) = /( j + | ) and //, (x) = /( |) to p/j = p;, = 8 bits, respectively 
(see Lemma 1). In the following bisection steps, we proceed in exactly the same manner. For instance, for the interval 
/ = ( |, j), we obtain j)(x) = {§g* 2 + — fjg which approximates /; to p; = 6 bits after the binary point. 

In a second step, we apply a modified Descartes method to isolate the roots of /. However, in- 
stead of computing the exact intermediate results obtained in the subdivision process, we only 
consider approximations to a certain number of bits. Whereas other methods [7,13,18,25,32] 
proceed in a similar way by using interval polynomials, our new method considers a specific 
approximation in each step and updates the possible approximation error. In [24,33], a similar 
approach was proposed. Therein, the proposed algorithms also initially start with an approxima- 
tion g of /, however, all intermediate results correspond to the initial approximation g and are 
computed exactly. In contrast, we propose to consider independent approximations of the inter- 
mediate results at each node of the recursion tree; see Figure 1.1 for a more detailed example. 

How is it possible that, for integer polynomials, an approximate version of the Descartes 
method is more efficient than the original "exact version"? Let us first consider the "exact 
Descartes method": Its complexity analysis shows that, for each interval (node) / = (a,b) in 
the recursion tree, the dominating costs are those for the computation of the Taylor expansion 
fi(x) := f(a + (b — a)x) at a; see Section 2.6 for a more comprehensive treatment. In each bi- 
section step, the polynomials /} ; and fj r (corresponding to the left and the right subinterval of /) 
are recursively computed from // by replacing x by x/2, followed by a Taylor shift by 1, that is, 
.14.1+1. More precisely, we have // ( (x) = //(x/2) and fi r (x) — fi t (x + 1 ) . In each iteration, the 
bitsize of the coefficients of // increases by n bits, and since the recursion tree has depth bounded 
by /Jmax = 0{wc), the representation of // eventually demands for at most T + nh max = 0{n 2 x) 
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bits. Hence, assuming asymptotically fast Taylor shift [15,40], the computation of a certain // 
amounts for (5(« 3 t) bit operations. 

Now, let us turn to the approximate method: In Section 2.3, we show that, for an arbitrary ap- 
proximation g of / to p max = d(nx) bits after the binary point, corresponding roots of / and g are 
almost at the same location with respect to their separations; see Theorem 3 and Appendix 6.2 for 
a more precise result. Thus, for each interval /, it should suffice to consider approximations // of 
// to p m ax bits after the binary point. Starting with an approximation of / to p max + 2h max = 0{m) 
bits after the binary point, we can iteratively obtain such approximations //. Namely, // can be 
recursively computed such that the approximation error quadruples at most in each bisection step 
and the height of the recursion tree is bounded by h max . Eventually, all polynomials // are rep- 
resented by 0{nx) bits (instead of 0(n 2 x) bits for the exact counterpart //) and, thus, the cost at 
each node decreases by a factor n. 

We will prove the above result for the more general setting where F is a polynomial with 
arbitrary real coefficients. More precisely, we show that it suffices to approximate each // to a 
number of bits after the binary point bounded by 0(L/ +n) = 0{Lf + nTp). Then, each fj is 
represented by 0(£f + T + nTp) bits and, as a consequence, the cost at each node is bounded by 
0(n(L F + T + nTp)) bit operations. We remark that, due to Appendix 2.3, we have T = 0(nTp) 
and, thus, the latter bound writes as 0(n{Lp + nYp)). The additional factor Lp + nFp in the 
bound (1.2) on the bit complexity is due to the size of the induced recursion tree. 

1.2. Outline 

In Section 2, we first introduce some basic notations. Furthermore, we derive a bound on how 
good / has to be approximated such that its roots stay at almost the same place with respect to 
the corresponding separations. Eventually, we revise the Descartes method before presenting our 
slight modification Dcm of it in Section 3. In Section 4, we present our new algorithm to isolate 
the roots of F and provide the corresponding complexity bounds. We conclude in Section 5. 
Parts of the complexity analysis as well as pseudo-code for our subroutines is outsourced to the 
Appendix. 



2. Preliminaries 

2.1. Some Notations 

For an interval / = (a,b), w(I) :=b — a denotes the width, m(I) := ^ the center, and r(I) 
" ! the radius of /. Furthermore, 

/ 4- ,+n / w(I) , , r , 7n / w(I) , w(I). 

denote extensions of / by and (to both sides), respectively. We will need these intervals 
for our modified version of the Descartes method as presented in Section 3. An (open) disc in 
C is denoted by A = A r (m), where meC indicates the center of A and r e R + its radius. The 
closure of a disc A or an interval / is denoted by A and /, respectively. 
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2.2. Scaling the Polynomial 



Instead of isolating the roots of the given polynomial F as in (1.1), we consider the equivalent 
task of isolating the roots of a "scaled" polynomial / which is defined as follows: We first com- 
pute an integer approximation T e N of the exact logarithmic root bound I> = log(max,' of 
F such that 

i> <r<4iog«+r>. (2.1) 

This computation can be done with O(inTf) 2 ) bit operations and demands for an approximation 
of F to 0{nYf) bits after the binary point; see Appendix 6.1. We can further assume that Y < T+l 
due to Cauchy's Bound [41] B C b '■= 1 +max, j^j < 1 + 2 T on the modulus of all roots. Now, we 
define 

f(x) = t^-= F{2 T' X) - (2-2) 

It follows that all roots z\ = %\ ■ 2~( r_1 \ . . . ,z n = %n ■ 2~( r_1 ) of / are contained within the disc 
Axrt(O) and the absolute value of each coefficient a,- of / is bounded by 2"( r+1 ' +T = 2°( nr \ In 
practice, it might be worth to investigate in an even tighter root bound F as described in [12, 
Section 2.4] in order to prevent the coefficients of / to become unnecessarily large. We further 
remark that the separations of corresponding roots of F and / scale by 2 r+1 (i.e., a(^,-,F) = 
2 r+1 -CJ(z 1 J)).Thus, 

E/ = -£loga(zi,/)=E F +n(r+l) = 0(nT + H F ). (2.3) 

i=\ 

2.3. Approximating Polynomials 

We assume that the coefficients of F are given as infinite bitstreams, that is, for a given pgN, 
we can ask for an approximation of F to p bits after the binary point. More precisely, each 
coefficient A,- is approximated by a binary fraction A,- = m,- -2~ p withm,eZand |A,-— A,-| <7r p ', 
e.g., A; = sign(A;)|jA;2''|j2~''. We call a polynomial F G Q[x] obtained in this way a p-binary 
approximation ofF. We remark that, in order to get a p-binary approximation of /, it suffices to 
approximate Fton(r+l)+p + T+ l bits after the binary point. Namely, given approximations 
A t =At + ^ with < pt := 2-("( r + 1 )+P+ T + 1 ) for all i, it follows that 



A L _A 1 

A n A n 



Mn+Mj < ^ . \ A n\ + \M < +2 r } < 2 -(„(r + l) + p) > 

|A n A n | \A n A n | 



Thus, a,- := j i -(2r)' approximates a,- = ^-(2 r+1 )' to an error less than 2 p . 

For an arbitrary polynomial g(x):= T.ILo gi x ' G C [x] with complex coefficients and an arbitrary 
non-negative real number fx e Rq , we define 

Mm := jlW = El^' e C[x] : |ft - 1,| < n for all i = 0, ... ,n| 

f/ze set of all ^-approximations of g. We remark that, since the coefficients of modulus less than 
jU can be approximated by zero, a jj. -approximation g of g might have lower degree than g. 
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Example. For g(x) := gfff^x 10 - 2x 2 + ^x - yg, the polynomial g(x) := ^x w - 2x 2 - jg con- 
stitutes a 6-binary approximation and g(x) := — 2x 2 — | a 2-binary approximation of g. 

2.4. Taylor Shifts 

For an arbitrary polynomial § e C[i] and arbitrary values m e C, A e M\{0}, let 

SKA]00 :=g(m + Ax). (2.4) 

The following lemma provides error bounds on how the absolute approximation error jj. of a 
polynomial g £ {g]^ scales under the transformation x i->- m + Ax: 

Lemma 1. For jU € Rq one/ g € [g]^ an arbitrary pt-approximation of a polynomial g £ C[i] of 
degree n, it holds that 

(in) ![_!,[] G k[_i i i ] ]2»/ 1 , anrf|[i,i] G fe[Li]]2" M - 

Proo/ For «(x) := (g — g) (x) = jU„x" + . . . + /iix + jUq, the absolute value of each coefficient 
is bounded by ju. Let m e C and A <G K\{0} be arbitrary values, then 

h(m + Xx) = £ W (in + Ax)' = £ W £ AV 4 ( ' ) = £ x* £ f ') (2.5) 

!=0 i=0 *:=0 Vv t=0 i=ifc vv 

Thus, for |m| < 1, the absolute value of the coefficient of x k is bounded by 

"W'-EH-'Q ^W'-EM''^) -mW-jt^, (2.6) 

where we used 

(i - M)-^ - E (" ( * + 1} ) (-i)'H' = E ("tO K - E (T) M'- 

!>o \ 1 / (>o V 1 / i>o V K / 

For m = A = 1/2, it follows that all coefficients of n are bounded by 2fX. This shows (i). For 
m = — an d A = 1 + jL, (2.6) implies that 

because * • ( l^}/^] ) < |j • sje < 4. Hence, (ii) follows. The first part of (iii) is also a direct 
implication of (2.6). The second claim in (iii) follows from the computation in (2.5) since each /x, 
is then (m = X = 1) bounded by „ • Q = L"=k (,_*) = EJS? ) < E£o* © < 2" ■ M- □ 

2.5. On Sufficiently Good Approximation 

In the next step, we derive a bound on how good / has to be approximated by an / such that, 
for all i, the distance of corresponding roots zi and z, of / and / is small with respect to the 
separation a(z,-,/). The following considerations are mainly adopted from our studies in [33]. 
Only for the sake of comprehensibility, we decided to integrate the results in this paper as well. 
We start with the following definition: 
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Definition 2. Let t > 1 be an arbitrary real value and / a polynomial as in (2.2). We define 



1 

- • min 

t i=\,...,n 



°(Zi,f)f'(Zi) 



8n 2 



We call a p 6 N sufficiently largJ^with respect to f if 

P>Pf= r-logAi(/,64« 2 )] = O^ + log/i-logKI) = 0(Z F +log«). 
The upper bound for P/ in (2.8) follows from 

°(zi,f)-\f(zi)\ = o(zuf)-\a n \H\zi-Zj\ > oizufi-WnlHoizjJ) = \a n \2- L f. 



(2.7) 



(2.8) 



and L f - log \a n \ = Z F +n(T+ 1) - log(2"( r+1 ) 
to our question raised above: 



= T.p. The following theorem gives an answer 



Theorem 3. Let f be the polynomial as defined in (2.2), t > 1 and f £ [/] 



Mz,J)/(tn) 



( i) For all i = 1 , . . . , n, the disc A; 
counterpart Zi of f. 

(ii) For each z 6 C\ U"=i A i> it holds that \ f(z) | > (n + l)n(f,t) 



contains the root Zi of f and a corresponding 



Proof. Since all roots of / are contained within Aj ^(O), it follows that a(zi,f) < 1 for all i and, 
thus, each disc A, is completely contained within the unit disc. For an arbitrary point z G d A; on 
the boundary of A, , we have 



l/WI 



> 



kini z ~ z /i : 

tn 

<y{z,J)\f{zi)\ 

tn 




II l*-*yl 

l<j<n.jjti 



°{z u f)\f{zi)\ 



1 

tn 



n-l 



> 



tn 

°(Zi,f)\f(Zi)\ 

2.12 -tn 



n 

l<j<nj^i 



■Zi 



■Zi 



>(n+l)n(f,t). 



In addition, since f e and \z\ < 1, we have |(/-/)(z)| < (w + l)ju(/,f) < |/(z)|. Hence, 

(i) follows from Rouche's Theorem applied to the discs A, and the functions / and /. For (ii), we 
remark that / is a holomorphic function on C\ U" =1 A; and, thus, \ f(z) \ becomes minimal for a 
point z on the boundary of one of the discs A, . □ 

From the last theorem, it follows that, for given / as in (2.2), it suffices to approximate the 
coefficients of / to p = 0{Lf + log« — log \a n \) = 0{Lf + logn) bits after the binary point to 
guarantee that each approximation / 6 [f]i-p nas its roots at almost the same location as /. 

Corollary 4. Let f be a polynomial as defined in (2.2) and p € N be sufficiently large with 
respect to f, that is, p > p/ with Pf as defined in (2.8). Then, each root Zi moves by at most 
when passing from f to an arbitrary approximation f G L/^-p- I n particular, real roots 
of f stay real and non-real roots stay non-real. Furthermore, for any z G C with \z — Zi\ > ^fir 
for all i, it holds that |/(z) | > (n + 1 )2~ p f . 



2 This definition is motivated by our results in Theorem 3 and Section 4. 1 
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2.6. The Descartes Method 



We first resume some basic facts about the Descartes method for isolating the real roots of 
a polynomial f(x) = Y!i=o a i x " <= Descartes' Rule of Signs states that the number var(/) 
of sign changes in the coefficient sequence of /, that is, the number of pairs (z, j) with ; < j, 
aiaj < 0, and a i+ i = . . . = ay_i = 0, is not smaller than and of the same parity as the number of 
positive real roots of /. If var(/) = 0, then / has no positive real root, and if var(/) = 1, / has 
exactly one positive real root. The rule easily extends to an arbitrary open interval / = (a,b) via 
a suitable coordinate transformation: The mapping x^ a + (b — a)x maps (0, 1) bijectively onto 
7, that is, the roots of / in / exactly correspond to those of 

//(*) := W)]W =f(a + w(I)x)=f(a+(b-a)x) (2.9) 

in (0, 1). Hence, the composition of x i->- a + (b — a)x and i^l /(l +x) constitutes a bijective 
map from (0, °°) to /. It follows that the positive real roots of 

fi^x) := (1 +*) n //(^y) = (1 +*)" -/(^) 

correspond bijectively to the real roots of / in 7. The factor (1 +x) n in the definition of fj lev clears 
denominators and guarantees that // rev is a polynomial. /}. rev is computed from // by reversing 
the coefficients followed by a Taylor shift by 1. We now define var(/,7) as var(// irev ). 

Based on Descartes' Rule of Sign, Vincent, Collins and Akritas introduced a bisection al- 
gorithm denoted Vca for isolating the roots of / in an interval 7o (here, we assume that 7o = 
(—1/2, 1 /2)). We refer the reader to [2,3,4,5,8,12] for extensive treatments and references. 



VCA. The algorithm requires that the real roots of / in Iq are simple, otherwise it diverges. In 
each step, a set srf of active intervals is maintained. Initially, srf contains 7q, and we stop as soon 
as A is empty. In each iteration, some interval 7 e srf is processed; If var(/,7) = 0, then 7 contains 
no root of / and we discard /. If var(/,7) = 1, then 7 contains exactly one root of / and hence is 
an isolating interval for it. We add 7 to a list & of isolating intervals. If there is more than one sign 
change, we divide 7 at its midpoint m(I) and add the subintervals to the set of active intervals. If 
m(7) is a root of /, we add the trivial interval [m(7),m(7)] to the list of isolating intervals. 



Correctness of the algorithm is obvious. Termination and complexity analysis of the VCA 
algorithm rest on the following theorem: 

Theorem 5 ([26,29]). Consider a polynomial f e M.[x], an interval I = (a,b) and v — var(/,7). 

(i) (One-Circle Theorem) If the open disc bounded by the circle centered at m(I) and passing 
through the endpoints of I contains no root of f(x), then v = 0. 

(ii) (Two-Circle Theorem) If the union of the open discs bounded by the two circles centered at 
m(I) ± z'(l/ (2\/3))w(7) and passing through the endpoints of I contains exactly one root of 
f(x), then v = 1. 

Proofs of the one- and two-circle theorems can be found in [2,12,22,26,27,28,29]. Theorem 5 
implies that no interval 7 of length oy or less is split. Such an interval, recall that it is open, 
cannot contain two real roots and its two-circle region cannot contain any nonreal root. Thus, 
var(/,7) < 1 by Theorem 5. We conclude that the depth of the recursion tree is bounded by 
1 / Of. Furthermore, it holds (see [12, Corollary 2.27] for a simple self-contained proof): 
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Theorem 6. Let I be an interval and h and h be two disjoint subintervals of I. Then, 

var(/,/ 1 )+var(/,/ 2 )<var(/,7). 

According to the above theorem, there cannot be more than n /2 intervals / with var(/, 7) > 2 at 
any level of the recursion. Therefore, the size of the recursion tree 7Vca is bounded by — nlog Of. 
For polynomials with integer coefficients of maximal bitsize T, it is shown that — logoy = 
0(n(\ogn + t)), thus, the latter bound writes as 0(n 2 x). However, a more refined argumenta- 
tion [12] shows that |7Vca| is even bounded by 0(nx). 

The computation of /}, rev at each node of the tree is costly. It is better to store with every 
interval I = (a 7 b) the polynomial fi(x) — f(a + x(b — a)). If 7 is split at its midpoint ra(7) into 
l[ = (a,m(I)) and I r = (m(I),b), the polynomials associated with the subintervals are fij(x) = 
//(f) and /,,(*) - f,C-f) = fi,(l +x). Also, / />rev (x) = (1 +x)"f I { T ^- x ). If the coefficients of 
/ are integers (or dyadic fractions) of bitsize T, then the coefficients grow by n bits in every 
bisection step. Thus, for a node 7 of depth h, the bitsize T/, of the coefficients of // is given 
bounded by T/j = T + nh. Hence, using asymptotically fast Taylor shift (see [40,15]), the number 
of bit operations needed to compute //,, fj r and //. rev from // is in 0(n(nh + t)). Since the depth 
of the recursion tree is bounded by 6{m), each // has coefficients of bitsize 0(n 2 x) and, thus, the 
cost at each node is in <9(n 3 r). Eventually, the total cost for VCA is in (5(h 3 t) • 6{nx) = 0{n A x 2 ). 

3. A Modified Descartes Method 

In the REAL-RAM model, where exact operations on real numbers are assumed to be avail- 
able at unit costs, the Descartes method can directly be used to isolate the real roots of the 
polynomial / as defined in (2.2). Then, for each node / of the recursion tree, we have to compute 
the number var(/,7) = var(// >rev ) of sign variations for the polynomial //, r ev an d the sign of / 
at the midpoint m(I). However, for an actual implementation, these computations turn out to be 
hard in general because the coefficients of / are arbitrary real numbers. To overcome this issue, 
we aim to only consider approximations of fi and //. rev instead. In Section 4, we will show that, 
for sufficiently good approximations of /, this approach is feasible. However, our approach does 
not directly apply to the Descartes method but to a slight modification of it. 

For our modified version of the Descartes method, we aim to replace the inclusion predicate 
var(/,7) = 1 by a predicate used in the Bolzano method; see Corollary 8. Section 3.1 resumes 
some useful results which are adopted from our studies on the Bolzano method [34] whereas, in 
Section 3.2, our modified version is formulated. 

3.1. The ,5^*(m, r)-Test: Existence of Roots 

For g e C[x], m e C and positive real values K and r, we consider the test 

,^(m,r): t s K {m,r) :=\g(m)\- 

k>l 

In order to simplify notation, we also write 3^* (A) or =5^(7) instead of ,5^*(m, r), where 
A = A, (m) or 7= (a,b) an interval with midpoint m = m(I) and radius r=r(I). If the polynomial 
g is fixed and no mix-up is possible, we further omit the "g" and write 5^(m,r) for 3?£(m,r) 

and 5^(m, r) for (m, r). We mainly use K = 3/2. Therefore, whenever the "K" is suppressed 
(i.e., we write 3? g (m,r) instead of £^y 2 (m,r)), we consider K = 3/2. Before presenting the main 

technical lemmata, we first summarize the following useful properties of 5^'(m,r): 



,(*) 



(m) 



k\ 



r k >0. 



(3.1) 
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• If &g(m, r) holds, then ST*, (m, r) holds for all K' < K and all r 1 < r. 

• For arbitrary values m, r and A ^0, the test ^(m,r) is equivalent to £? K ' m ' A ' (0,r/A) be- 
cause of f jr [m '^ 1 (0,r/A) = ?£(m,r). In particular, for an interval / = (a,fe), the test ^'(0,r) is 
equivalent to £?£(a,rw(I)), where g/(x) = g(a + w(/);c). 

• For A e M + , ff-(m,r) = t^ 8 (m,r) ■ A -1 and, thus, £^(m,r) is equivalent to &£ g (m,r). Hence, 

''(m,r) and (m,r) are equivalent since (#/)' = (g(a + w(I)x))' = w(I)(g')i. 

The 5^?(m,r)-test serves as exclusion predicate but might also guarantee that a certain disc 
contains at most one root. We refer to [6, Theorem 3.2] for a proof of the following lemma. 

Lemma 7. Consider a disc A = A m (r) C C and a polynomial g G R[x\: 
(i) If £?k{A) holds for a K > 1, //zen A contains no root of g and 

(1-I)| g („)|<| g ( z )|<(l + I)| g („)| 

/or z z'n f/ze closure A o/A. 
(») If ^fy 2 (A) holds, then A contains at most one root of g. 

The 3^y 2 (m, r)-test now easily applies as an inclusion predicate: 

Corollary 8. Let I — (a, b) be an interval such that 5^' 2 (0,r) holds for an r> 1. Then, I contains 
a root t; of g exactly if g(a) ■ g(b) < 0. In the latter case, the disc A rw ^) (a) is isolating for t;. 

Proof. If &y' 2 (0,r) holds, then ^/ /2 (a,rw(I)) holds as well. It follows that the disc A rH , (/) (a) 
and, thus, / contains no root of the derivative g' . Now, since / is monotone on /, it suffices to 
check for a sign change of g at the endpoints of /. Namely, there exists a root | of g in / if and 
only if g(a)g(b) < 0. In case of existence, A w (,j (a) is isolating for E, due to Lemma 7. □ 



In order to show that the iJ^ 2 (m, r)-test in combination with sign evaluation is an efficient 
inclusion predicate, we give lower bounds on r in terms of a g such that the predicate succeeds 
under guarantee. 

Lemma 9. For g a polynomial of degree n, a disc A = A r (m) C C, an interval I — {a,b) and 

I+ = ( fl - ^ir > b + )' u holds that: 

(0 Ifr < 0, then ST(A) or ST' (A) holds. 

(ii) If A contains a root £, of g and r < then i^'(A) holds. 

(Hi) Ifvax(g,I + ) > and £? g 'i(0,2) fails, A 2w (/)(a) contains a root E, of g with a(B,,g) < 8n 2 w(/). 
(iv) Ifvar(g',I) > and £? S '(Q, 1) fails, A 2nw ^(a) contains a root t, of g with o(E,,g) < 4n 2 w(I). 

Proof. For the proof of (i) and (ii), we refer to [33, Lemma 5]. For (hi), suppose that var(g,/ + ) > 
and £T g i(0,2) does not hold. Then, according to Theorem 5 (i), the disc A w ^+y 2 (m(I)) C 

A 2H ,.(/)(a) contains a root £, of g. With (ii), it follows that 2w(I) > a ^ n i^ and, thus, <j(E,,g) < 
8n 2 w(/). For (iv), we first argue by contradiction that A2„ w (/)(a) contains a root ^ of g: If \a — 
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Xi\ > 2nw(I) for all roots x,- of g, then 

r 1 < (r — V < (— V 

where the prime means that the ;' 7 's (J=l...k) are chosen to be distinct. It follows that T g (a,w(Ij) 







8(a) 





holds because of YH=\ ^jn?r w (I) k < YX=\ ^ * < 1 < §• In addition, Theorem 5 guarantees the 
existence of a root EJ <E A^mnfrnf/)) of g'. Hence, we have |<i; — <i;'| < 2nw(I) + r(I) < 4nw(I) 
which implies u(^,g) < 4n 2 w(I) due to the fact [11,41] that there exists no root of the derivative 

s'inA CT( ^ ))Al (£). □ 



3.2. DCM: A Modified Descartes Algorithm 

We introduce our modified Descartes method DCM (short for "Descartes modified") to isolate 
the real roots of a polynomial / as defined in (2.2). We formulate the algorithm in the REAL- 
RAM model, thus, it still does not directly apply to bitstream polynomials. However, in Sec- 
tion 4.1, we will present a corresponding version DCM P of DCM which resolves this issue; see 
also Appendix, Algorithm 1 for pseudo-code of DCM. 



DCM. DCM maintains a list srf of active nodes and a list of isolating intervals, where we 
initially set — and srf := {(/o,// )} with := (— \, j). For each active node (/,//) £ 
we proceed as follows. We remove (/,//) from si '. Then, we compute the number Vj+ := 
var(/,/ + ) = var(/ 7 + rev ) of sign variations for / on the extended interval /+. We remark that 
//+ to = fi(-l + (1 + sW and / /+>rev (x) = (1 +x) n f I+ (^- x ). If v /+ = 0, we do nothing. If 
V/+ > 1, we consider the test 3Th (0,2) which is equivalent to (a,2w(I)). If it fails, then / is 
subdivided into // = (a,m(I)) and I r = (m(I),b) and we add (//,//,) = (Ii,fi(%)) and (I r ,fi r ) = 
(I r ,f Il (x+ 1)) to srf. Otherwise, we evaluate the sign s of f(a + ) ■f [b + ) = f I+ (0) If s < 

and I + is disjoint from any other interval in G, we add /+ to 0. If s > or / intersects an interval 
in G, we do nothing. The algorithm stops when srf becomes empty. 



Theorem 10. For the polynomial f as defined in (2.2), DCM terminates and returns a list = 
{h,. ■ ■ ,I m } of disjoint isolating intervals for all real roots off. 

Proof. If the width w(I) of an interval / = (a,b) is smaller or equal to ^4, then, according to 

Theorem 5, var(/,/ + ) = or 3^^(0,2) holds. Thus, / is not further subdivided. This shows 
termination of DCM. From our construction and Corollary 8, each interval in is isolating for 
a real root of / and all intervals in are pairwise disjoint. It remains to show that, for each real 
root E, of /, there exists a corresponding isolating interval in 0. Since all roots of / have absolute 
value bounded by 1 /2, there must be a terminal interval / = (a,b) whose closure I contains %. 
Since Vj+ > 0, / cannot be discarded in the first step of DCM. Hence, !?fi (0,2) holds and, thus, 
/ is monotone on / + . Since / + contains the root E,, we have f(a + )f(b + ) < 0. It follows that 
either /+ is added to the list of isolating intervals or I + intersects an interval J + = (c + ,d + ) € 
which has been added to before. Let J = (c,d) be the corresponding smaller interval for J + . 
Since the ^^-neighborhood of / intersects the ^^-neighborhood of J, the preceding Lemma 1 1 
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shows that one of the discs A 2vv (/)(a) or A 2vl ,(/)(c) contains both intervals 7 + and J + . Since both 
i^r'(0,2) and ^^(0,2) hold, each of the latter two discs contains at most one root due to 
Corollary 8. It follows that J + G & already isolates ^. □ 

Lemma 11. [^Lef / = (a,b) and J — (c,d) be two intervals (not necessarily of equal length) 
of the form (-§ + i2~*,-| + (i+ l)2- h ) , where h G N and i G {0,...,2 h - 1}. If the 
neighborhood £/ w (/)/2nC0 °fl intersects the "^--neighborhood U w fj\n n (J) of J, then one of the 
discs A 2w (/)(<z) or ^2w(J) ( c ) contains the intervals (a — w(I),b + w(I)) and (c — w(J),d + w(J)). 

Proof. W.l.o.g., we can assume that w(J) > w(I) and, thus, w(J) = 2 l w(I) with an / € No- Let 5 
denote the distance between 7 and /. If 8 = 0, then A 2w ( /) (c) contains (a — w(I) ,b + w(I)) and 
(c - w(J) , d + w(J) ) . If 8 ^ 0, then 8 = 2 k w(I) with ateN . Since U w{1) /2 „ (/) n U w(J) /2 „ (/) ^ 0, 
we must have > f . In particular, we have > | = 2 k ~ l w(I). Since w(7) and w(J) differ 

by a power of 2, it follows that w(J) > 2 k+2 w(I) = 48 and, thus, 2w(J) = w(J) + ^ + ^ > 
w{J) + 2w(I) + 28. From the latter inequality our claim follows. □ 

Theorem 12. For a polynomial f as in (2.2), Dcm induces a subdivision tree 7d C m of 

height /z(Tdcm) = O (log n — log of) and size |7dcm| = 0(E/ + nlogn). 

Proof. The result on the height of Tqcm follows directly from the proof of Theorem 10. Namely, 
we have shown that Dcm never subdivides an interval of width less than or equal to r^. For 
the bound on |7dcm|> we use a similar argument as in [14] and [24]. Namely, for a root £ of 
/ and a certain h G N we say that I = (-± +/2 _/, ,-j + (; + l)2" /! ), / = {0, . . . ,2' ! - 1}, is a 
canonical interval for % if the real part of ^ is contained in [— i +i2~ h , — j + (; + l)2 _/l ) and 
cr(|,/) < 8n 2 2 _,! = 8« 2 w(7). We denote T c the canonical tree which consists of all canonical 
intervals. We remark that, for a canonical interval 7, the parent interval of 7 is canonical as well. 
The following considerations will show that |7dcm| = 0(|T e |) and \T C \ = 0{Lf + «log«). For the 
size of the canonical tree, consider a leaf 7 G T c and let ^/ be a root of / corresponding to this 
leaf. If there are several, then is the root with minimal separation. Then, a(§/,/) < 8n 2 2~ /! 
and, thus, h < 21ogn + 4 — log a(§/,/). Since each root of / is associated with at most one 
leaf of the canonical tree, we conclude \T C \ = 0(n\ogn + £/). It remains to show that |7dcm| = 
0(|T C |). Consider the following mapping of internal nodes (intervals) of Tqcm to canonical nodes 
(intervals) in T c : Let 7 be a non-terminal interval of width w(I) = 2~ h . Then, var(/,7 + ) > 
and T-fi (0,2) does not hold. According Lemma 9 (iii), the disc A 2vv (/)(a) contains a root t, of / 
with ct(§,/) < 8n 2 w(7) = 8n 2 2~ /! . Hence, one of the four intervals l\ = (a — 2w(I),a — w(7)), 
7 2 = (a — w(I), a), I or 7 2 = (b, b + (b — a)) is canonical for | . We map 7 to the corresponding 
interval. This defines a mapping from the internal nodes of 7b CM to the nodes of the canonical 
tree T c . Furthermore, each node in the canonical tree has at most four preimages in Tdcm and, 
thus, the number of internal nodes of Tqcm is bounded by 0(n\ogn +£/). Since Tdcm is a binary 
tree, the bound on the number of internal nodes applies to the whole tree as well. □ 



3 Lemma 1 1 proves a slightly stronger result than necessary for the proof of Theorem 10. The stronger result applies in 
the proof of Theorem 15 in Section 4.2. 
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4. Algorithm 

We first outline our algorithm RISOLATE to isolate the roots of /. BTSOLATE decomposes into 
two subroutines Dcm p and Certify'', where p indicates the actual working precision. Dcm^ is 
essentially identical to Dcm with the main difference that, at each node / = (a, b) of the recursion 
tree, we only consider approximations fi(x) of fi(x) = f(a + w(I)x) to a certain number p/ of 
bits after the binary point, where p + 21ogw(/) < pi < p. We remark that we proceed / in a way 
such that it is terminal for Dcm'' if it is terminal for the exact counterpart Dcm. This ensures 
that, for any p, Dcm'' induces a subtree T Dcm p of Tucm and, thus, |r DCM p| = 0{Lf + nlogn) 
due to Theorem 12. We further show that, for a precision p > p™ ax = 0(Lf + n), Dcm'' returns 
isolating intervals for all real roots of /; see Theorem 15 for the definition of p™ ax and further 
details. However, for smaller p, Dcm'' may return isolating intervals only for some roots but 
without any information whether all real roots are captured or not. In order to overcome such an 
undesirable situation, we consider an additional subdivision method Certify'' similar to Dcm'' 
which aims to certify that all roots are captured. We further show that CERTIFY* 5 also induces a 
recursion tree of size 0(E/ + nlogn) and succeeds if p > p™ ax . If, for a given precision p, our 
algorithm fails to isolate all roots of /, we double p and restart. 

4.1. Dcm'': An Approximate Version o/DCM 

We present our first subroutine Dcm''. Comments to support the approach are in italic and 
marked by a "//" at the beginning. 



DCM P . Let Iq = (— j) 3) b e me startm g interval which, by construction of /, contains all real 
roots of /. In a first step, we choose a (p + n + 1) -binary approximation / of / and evaluate 
/(—j +x). Then, the resulting polynomial is approximated by a (p + 1) -binary approximation 

//„ G [/(— 5 +x)] 2 -p-i and, according to Lemma 1, we have fj <E \fi ]2-p- 

DCM P maintains a list s/ of active nodes (/,//, pi), where / = (a,b) C Iq is an interval, // 
approximates // to p/ bits after the binary point and p +21ogw(7) < Pi < p. DCM P eventually 
returns a list of tuples (J,sji,sj r ,Bj), where J = (c,d) is an isolating interval for a root 
°f /> s J,i — sign/(c), sj tr = sign/(J) and < Bj < min(|/(c)|, \f(d)\). We initially start with 
si := {(I(),fi ,L)} and & := 0. For each active node, we proceed as follows: 

(1) Remove (/,//, p/) from si . 

(2) Compute the polynomials 

Mx)=M~ + (l + ^ i )x) and h(x) = ph i x i -.= (l+x) n f I+ (j^). (4.1) 

(3) If k > -2 n+2 -P' for all i or n, < 2 n+2 -P> for all i, do nothing (i.e., / is dicarded). 

II A simple computation (see the subsequent Lemma 13 (i)) shows that h approximates 
//+,rev( x ) = (x+ l)"//+( ]4^) to Pi — n — 2 bits after the binary point. Thus, i/'var(/,/ + ) = 
0, all coefficients ofh are either smaller than 2 n+2 ~ pl or larger than —2 n+2 ~ p] . Since we 
want to induce a subtree of the recursion tree Tq cm induced by f, we discard I if all coef- 
ficients ofh are larger than —2" +2 ~P' (or smaller than 2 n+2 ~ Pl ). 
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(4) If there exist k and hj with k < -2 n+2 -P' and hj > 2 n+2 -P', consider the test ST^' (0, 2), 
that is, evaluate t^ (0,2). 

// Due to Lemma 13 (i), we have \t$ (0,2) -?^'(0,2)| < n2 n+l -P'. Hence, if (0,2) 



'3/2 '3/2 

holds, then ty^ (0,2) > — n2 n+l ~ p <. Thus, we proceed as follows: 



(a) If t^'(0,2) > -n2 n+l -P', consider the polynomial 

f I (x):^f I (x)+n2 n+l -P'-x, (4.2) 
// Then, Sf^M (0,2) holds and, in particular, fj is monotone on (—2,2). 



evaluate 



A-:=//(-^)=/ /+ (0)-2"- 1 -"' (4.3) 
A+ := f,(\ + 1) = / /+ (1) + (An + l)2"- l -P> (4.4) 



4n 



A :=//(--)=//(--) -2" +1 -"', 



(4.5) 



and check whether the following conditions are fulfilled: 



w(I) w(I) 

1= (a,b) = (a- ——,/?+——) intersects no J for any (J ,sj i,sj^ r ,Bj) G 6, (4.6) 
In In 

A~-A + <0, (4.7) 

min(|A-|,|A+|) >2 n+3 - p 'n, and (4.8) 

|^| > 2 deg//+«+7-p, w 2_ (4 9) 

If any of the conditions (4.6)-(4.9) fails, do nothing. If all conditions are fulfilled, then 

add (7,signA-,signA + ,min(|A-|,|A+|)-2' 1+3 -'''n) to G. 

II If (4.7)-(4.9) hold, then I is isolating for a root of f; see Lemma 13 (Hi). Further- 
more, since fj is monotone on (—2,2), we have \fi(— ^)| > |A~| and |//(1 + ^)| > |A + |. 
Then, from inequality (4.8) and Lemma 13 (ii), it follows that sign / '(a) = sign (A ~), 
sign/(£)=sign(A+) andmm(\f(a)\, \f(b)\) > min(|A-|, |A+|) -2"+ 3 -<"«. 

(b) If t^jl (0,2) < -n2 n+l -P', subdivide / into /, := (a,mf) and I r := (m u b). Compute a 

P/-binary approximation // ( of //(f) and a (p/ — l)-binary approximation fi r of fi(^-), 
and add (//,// ; ,p/— 1) and (/ r ,// r ,p/ — 2) to .e/. If p/ < 2, return "insufficient precision". 

// Due to Lemma 1, we have fj t G L//J2-P/- 1 and f/ r G [//,.] 2 -P/- 2 - Hence, by induction, it 
follows that p + 2 log w(/) < p/ < p /or aZ/ acfi've nodes. 

Dcm'' stops when «e/ becomes empty. It may either return "insufficient precision" (in Step 4 (b)) 
or a list of isolating intervals / for some of the roots of / together with the signs of / and a 
lower bound on |/| at the endpoints of I. 
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A 2 (0) 



Fig. 4.1. If A~ ■ X + = -fi(l + <0, then there exists a root y e (- ^ , 1 + ^) of /;. Furthermore, A 2 (0) 

contains no further root of fj. A computation shows that |/;(z)| > |//(z) — //(z)| for all z on the boundary of the 
~ -neighborhood 1/ of (0, 1) if the inequality (4.9) holds. Then, due to Rouche's Theorem, U isolates a root of /;. 

Lemma 13. Let f be a polynomial as in (2.2), I = (a,b) an interval considered by DCM^ and h 
the polynomial as defined in (4.1). Then, 

(i) h{x) £ [f I+ , tev ] 2 n+2- Pl and \t$(0a)-4 f /}(0,2)\<n-2" +1 -P'. 

(ii) For an arbitrary real value t with \t \ < 1 + ^, it holds that \ f(a + 1 ■ w(I)) — fi(t) \ < 2 n+i ~P'n, 
with fi as defined in (4.2). In particular, 

|/(a+)-A-|, \f( a -^l)-X\, \f(b + )-X + \<2 n+i -P>n, 
n 

with A~, A + and A as defined in (4.3)-(4.5). 
(Hi) Suppose that ty'2 (0,2) > — n2 n+1 ~ p > and the inequalities (4.7)-(4.9) hold. Then, I + contains 

a real root t, of f and the ^^--neighborhood of I is isolating for 
(iv) For any tuple (J,sj j,sj , r ,Bj) € G, the endpoints of J are located outside the union of the discs 

Ai := A o-( z „/)/(64«3)fc)> where i=l,...,n. 

Proof. Since // £ [//]2-P/> we nave .//+ <= [//+12-P/+ 2 due to Lemma 1 (ii). Reversing the coeffi- 
cients and replacing x by x + 1 increases the error by a factor of at most 2" (see Lemma 1 (iii)), 
thus h £ [//+,revj2-P/+ 2 +"- F° r t ' le secon d P art of (i), consider the following simple computation: 

\t§[ (0, 2) - rg' (0,2)| < | -n ■ 2-P< £ 2 l = \-n- T* (2" - 1 ) < nT+^P' , 

where the first inequality uses (fi)' £ [(./})'] n -2~P/ ■ 
For (ii), we have 

i/(«+fw(/))-/ / (oi-i//(o-//(oi<i//W-//(oi+ki-2" +1 - p '» 

< 2-P> £ |»|' + (1 + -)2 n+l -P'n < n2-P'(\ + +n2" +2 -P< < n2 n+3 -P> 

i=Q n 

Now, if the inequalities (4.7) and (4.8) hold, then sign/(a + ) = sign(A~), sign/(£> + ) = sign(A + ) 
and f(a + ) ■ f(b + ) < 0, hence, / has a real root in I + . We next show that (4.9) implies the 

( f V ( f V 

uniqueness of this root. From t%U (0,2) > -n2 n+1 ~P', it follows that ^ /2 n (0,2) succeeds 
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and, thus, A2(0) contains at most one root of //. Since A = //(— i) and A + = //(l + 
have different signs, the interval (— ^,1 + ^) contains a root 7 of //. We consider the \- 
neighborhood U C C of (0, 1) and an arbitrary point z on its boundary; see Figure 4.1. It holds 
that I - \ -Ti/lz-yj < (1 + ^)/(^) =4n + 3 < 8n and, for any root 77^ 7 of //, we have 

l-^-y| < l-^-^l + k-rl <1 i + ^ =2 i + a 

Hence, it follows that 



//(z) 



li ll TT I " f\ < n 2 deg ^ /+2 (l + J_) de s//-l(l _ I)-deg//+l 



m 

< n2 de ^'+ 2 • V2?72 • 2.72 < «2 de s / '+ 4 



and, thus, |/,(z)| > |A| • 2- de s^- 4 n -\ Since |z| < 1 + \, we have |//(z) -/,(z)| < n2 n + 3 ~P' 
according to (ii). Then, from Rouche's Theorem, it follows that // has exactly one root within 
U if (4.9) holds. This shows (iii). It remains to prove (iv): Let J = 1 = (a,b) and / = (a,b) the 
corresponding smaller interval. From our construction and (iii), /+ contains a root E, = Zi of / 
and the '^-neighborhood of / is isolating for this root, thus, \a — Zi\ > for all i. If there 
exists an ;' 7^ z'o with a G A,-, then w(I) < 4n\d — Zi\ < <*(zi,f)/ (16n 2 ). Thus, we obtain 

l 5 - a |<l 5 - a | + | a - a |<( 1 + i w/) + ^<( 1 + i)£|^ + ^)< ofc , / ), 

a contradiction. It remains to show that a £ A, . If a e A, , then w(I) < ■ According to 

Lemma 9 (ii), T^' (0, 2) already holds for a parent node J of / and, thus, t^J (0, 2) > n2 n+l -P> 
because of (i). This contradicts the fact that J is not terminal. In completely analogous manner, 
one shows that b is also not contained in any A,. This proves (iv). □ 

We close this section with a result on the size of the recursion tree induced by Dcm^ and the 
bit complexity of Dcm p : 

Theorem 14. Let f be a polynomial as in (2.2) and p G N an arbitrary positive integer. Then, 
the recursion tree T Dcm p induced by Dcm p is a subtree of the tree Tq C m induced by DCM, thus, 

\ t dcmp\ < |7dcm| = 0(£/ + «logn). 

Furthermore, DCM P demands for a number of bit operations bounded by 

6 (n (L f + log n) (nT + t + p - log a f ) ) . 

Proof. For the first claim, we remark that Dcm'' never splits an interval / which is not split 
by Dcm when applied to the exact polynomial /. Namely, if / is terminal for Dcm, then ei- 
ther f|/2 (0; 2 ) > or var(/,7+) = var(/ /+ rev ) = 0. In the first case, we must have t^j\ (0,2) > 
— n2 n+l ~ Pl whereas, in the second case, all coefficients hi of h(x) = (1 + x) n fj+ (]4^) are either 
larger than — 2 n+2 ~P' or smaller than 2 n+2 ~ pl ; see Lemma 13 (i). Thus, / is terminal for Dcm^ 
as well. The result on the size of T Dcm p then follows directly from Theorem 12. 

For the bit complexity, we first consider the cost in each iteration: For an active node (/,//, pi) € 
srf, I = (a,b), the polynomial // approximates // to p/ < p bits after the binary point. The 
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absolute value of each coefficient of // is bounded by 2 n+x (2T) n because the shift operation 
x i y a ~\~ (b a)x does not increase the coefficients of / by a factor of more than 2" and the ab- 
solute value of the coefficients of / is bounded by 2 T+ "( r+1 ); see Section 2.2. It follows that the 
bitsize of the coefficients of // is bounded by n(T+l) + T + p. Hence, the cost for computing 
h(x), fi t and fi r (x) is bounded by 0(n(nT + x + p)). Namely the latter constitutes a bound on 
the cost for a fast asymptotic Taylor shift by an O (log n) -bit number. The cost for evaluating 
4/2 (0>2), A - , A + and A matches the same bound because all these computations are evalua- 
tions of a polynomial of bitsize 0(nY + x + p ) at an (9(log n)-bit number. We further remark that, 
in each iteration, 6 contains disjoint isolating intervals / for some of the real roots of / and, 
thus, \0\ <n. Hence, the endpoints of the interval J have to be compared with those of at most 
n intervals stored in . Since Dcm p does not produce any interval of size less than ^4, these 
comparisons demand for at most O (n (log n — log oy)) bit operations. It follows that the total cost 
at each node is bounded by d(n{nT+T + p — logo/)) bit operations. The bound on the total cost 
then follows from our result on the size of the recursion tree. □ 

4.2. Known Lf and Of 

From Corollary 4, we already know that, for p > p/, each root z, of / moves by at most a ^ n P 
when passing from / to an arbitrary approximation / £ [/] 2 -p/ '■> see Definition 2 for the definition 
of Pf. Hence, we expect it to be possible to isolate the roots of / by only considering approx- 
imations of / (and the intermediate results //) to Pf bits after the binary point. The following 
theorem proves a corresponding result. 

Theorem 15. Let f be a polynomial as in (2.2) and p G N an integer with 

p>pp x :=\p f -3\ogGf + 16n\ =0(E/ + n). (4.10) 
Then, DCM P returns isolating intervals for all roots offandBj > 2 p f for all (J,Sjj,Sj t r,Bj) £ G. 

Proof. Due to Theorem 12 and 14, the height /z(Dcm p ) of r DCM p is bounded by 

Ii(Dcm p ) < log^- =21ogn + 4-logoy <4n- logoy. 

Then, for any interval I = (a,b) produced by Dcm p , we have 

p / >p+21ogw(/)>p-2/!(DcM'')>p^ in := [p 7 + 8n- logoy] > 0. (4.11) 

The latter inequality guarantees that DCM* 5 does not return "insufficient precision". Now let / be 
an interval whose closure I contains a root E, = Zi of /. We aim to show the following facts: 

(1) / is not discarded in Step 3 of DCM P . 

(2) If 4^' (0,2) > -n2 n+l -P', then all inequalities (4.7)-(4.9) are fulfilled. 

(3) In the latter case, either I ={a— , b + ) is added to & or / only intersects intervals 
7, with a corresponding (J,sjj,sj t r,Bj) £ 6, which already isolate E,. 

If (l)-(3) hold, then DCM P outputs isolating intervals for all real roots of /. Namely, DCM P 
starts subdividing /q = (— j, 5) which contains all real roots of /. Thus, for each root B, of /, we 

— if v - 

eventually obtain an interval / such that / contains | and t^U (0,2) > — n2 n+ pl . Then, either / 

is added to the list of isolating intervals or & already contains an isolating interval for E, . 
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For the proof of (1), we have already shown that w(I) > ■ Corollary 4 then ensures that 
an arbitrary g € [/] 2 -p/ nas a root <■>' <= I + ■ Namely, the root | e 7 stays real and moves by at 
most CT i|^ < when passing from / to g. Now, suppose that all coefficients hi of ft (x) = 
(l+x)"/ /+ ( r L) are larger than -2" +2 -P'; see (4.1) for definitions. Since \hi~hi\ <2 n+2 -P' for 
all coefficients h, of //+ rev = L" =0 ft;*' (see Lemma 13 (i)), it follows that hi > —2 n+3 ~ p > for all 
i. Hence, for the polynomial 

g(x):=f(x)+2" + '-P'e[f} 2 - Pf , 

we have gi+ :KV (x) = fi+,rev( x ) + 2" +3 ~ L/ (x + 1)" and, thus, g J+ rev has only positive coefficients. 
In the case where hi < 2 n+2 ~ pl for all i, we consider g(x) := f(x) - 2" +3 ~ Pl e [/] 2 -p/ and, thus, 
g!+ p has only negative coefficients. Hence, in both cases, there exists age [/] 2 -p/ which has no 
root in /+, a contradiction. It follows that / cannot be discarded in Step 3. 

For (2), suppose that 4/ 2 ' (0,2) > -n2 n+1 -P'. Due to Lemma 13 (i), we have t^' (0,2) > 
-n2 n+2 ~ pl , and since log " 2 " w + ( ^ P/ < 6 + 3 log n + n - p/ - log oy < -p/, it follows that 

n2 n+2-p, 
:= /(*) +x- e [/] 2 - P/ . 

Hence, g has a root |' in /+. Since 4/2 '(°> 2 ) = 4/2 '(°' 2 ) + » 2 " +2 ~ L/ > 0, the disc A 2w(/) (a) 
is isolating for The following argument shows that A iw ^y 2 (a) isolates E,: Suppose that 
A 3w (/)/ 2 (fl) contains an additional root z ; - ^ <!; of /. Then, < 3w(7) and, thus, <^ and 

Zy would move by at most < when passing from / to g. It follows that g would have at 
least two roots within A 2w (/)(a), a contradiction. Now, since A 3w (/)/ 2 (a) is isolating for % € 7, we 
have 

° ( ^ f) <w(I)<2o($,f). 



16n 2 

The left inequality implies that the distance of | to any of the points a + =a — v ^-, b + = b + 



and c := a — is larger than or equal to > ■ Let di :— \zi — a\ denote the distance 
between a root Zi 7^ % and the disc A 3w ( / )/ 2 (a). Then, 

°(zj,f) < \zj-Z\ < 4 + 3w(7) w(7) 
64n 3 ~ 64n 3 ~ 64n 3 ' 4 

It follows that the points a + , b + , c € A5 W (/)/4(a) are located outside the disc A; := A^. /)/(64„3)(zi)- 
In summary, none of the discs A„ i = l,...,n, contains any of the points a + , b + and c. Hence, 
due to Corollary 4, it follows that each of the values |/(c)|, \f(a + )\ and \f(b + )\ is larger than 
(n+ l)2~ p f . A simple computation now shows that (n + l)2 _p / > 2 2 " +8_P/ ra 2 . Thus, according 
to Lemma 13 (ii), each of the absolute values |A|, |A~ | and |A + | is larger than 

(n + 1)2- p f -2 n+i - p 'n > 2 2n+% - p 'n 2 - 2 n+3 - p 'n > 2 2n+1 - p 'n 2 . (4.12) 

It follows that the inequalities (4.8) and (4.9) hold. Since I + is isolating for f(a + ) and f(b + ) 
must have different signs and, thus, the same holds for A~ and X + . Hence, the inequality (4.7) 
holds as well. In addition, we have Bj = min(|A _ |, |A+|) -2 n+3 - pI n > 2~ p f because of (4.12). It 

remains to show (3): If ty 2 (0,2) > —n2 n+ p ', then due to (2) and Lemma 13 (ii), the interval 
I and the ^^-neighborhood of I is isolating for If /does not intersect any other interval in 0, 
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then / is added to and, thus, DCM P outputs an isolating interval for % . We still have to consider 
the case where / intersects an interval J from 0. From the construction of ff, J is the extension 
(c, d) of an interval f = (c,d). Now, suppose that J is isolating for a root 77^^. The roots | and 7 
move by at most and > ^p, respectively, when passing from / to an arbitrary g G [/] 2 -p/ ( see 
the proof of ( 1 )) . Hence, it follows that the union of (a — w (I) , b + w(I) ) and (c — w(J'),d+w(f)) 
contains at least two roots of any g G [f] 2 -p f . Due to Lemma 11, one of the discs A 2w ^(a) or 

/„ y 

A 2m ,(//)(c) then also contains at least two roots of g contradicting the fact that ty' 2 (0,2) > for 

p(x) :=f( x )+ x .^ffL g [/] 2 _ p/ and ^'(0,2) >0far«(x) :=/(*)+*■ G [/] 2 - P/ . 

It follows that J already isolates E, . □ 

4.3. Unknown Pf and 07 

For unknown Pf and Of, we proceed as follows: We start with an initial precision p (e.g., 
p = 16) and run DCM P . If DCM P returns "insufficient precision", we double p and start over. 

Otherwise, Dcm p returns a list & — {(Jk,Sk,h s k,r,Bk)}k=i m, where each interval — (ck,dk) 

isolates a real root of /, s k j = sigaf(c k ), s k , r = sign/(4) and < B k < min(|/(c t )|, \f{d k )\). As 
already mentioned, there is no guarantee that all roots of / are captured. Hence, in a second step, 
we use the subsequently described method Certify p to check whether the region of uncertainty 



1 1 

"2'2 



k=i 



may contain a root of /. If we can guarantee that f(x) ^ for all x G we return the list 
Jz? = {Jk}k=\ m °f isolating intervals. Otherwise, we double p and start over the entire algo- 
rithm. We have already proven in Theorem 15 that Dcm'' isolates all real roots of / if p > p™ ax 
(i.e., p fulfills the inequality (4.10)). The following considerations will show that, for p > p™ ax , 
Certify^ succeeds as well. 

How can we guarantee that / does not vanish on Ml The crucial idea is to consider a decom- 
position of [—5, 5] mt0 subintervals / and corresponding ^/-approximations g of // such that g 
is monotone on [0, 1] or £^y 2 (0, 1) holds. Namely, for such an interval /, we can easily estimate 
the image g([0, 1]) and, thus, conclude that / contains no root in in M or p < pj^ because g{t) 
and //(f) differ by at most (n + l)jti/ for all t G [0, 1]. More precisely, we have: 

Lemma 16. Let I = (a,b) be an interval and g(x) a ^.-binary approximation of // with 

-logju >p-2(4n-logoy). (4.13) 

(i) Suppose that i£^ 2 (0, 1) holds and I is not entirely contained in one of the J k . If 

|^(0)| > 8«M, (4-14) 



then I contains no root of f. Otherwise, p < p™ ax . 
(ii) Suppose that g is monotone on [0, 1] and let ID 2% = |J;=i ^; be the intersection of I and . 
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Fig. 4.2. DcmP returns a list G = {Jk,>>kji s k.r-Bk}k> where is isolating for a real root of /, s^j = sign/(cj.), 
s k.r = sign/t^t) an d m i n (l/( c /t)l > l/(<4)l) > > 0. The intervals in between define the region of uncertainty 2%. In 
CERTIFY'', we subdivide ( — 1/2, 1 /2) into intervals I such that, for a ^-approximation g of //, either 5j^ 2 (0, 1) holds 
or g is monotone on (0, 1). If 5^ 2 (0, 1) holds and |g(0)| > Smfl, then / contains no root of /; see Lemma 16 (i). If g 
is monotone on (0, 1), we consider all intervals L, in the intersection of / with M and check whether the conditions in 
Lemma 16 (ii) are fulfilled. If they are fulfilled, then / has no root in L, ; otherwise, we must have o~ < a'™ x . 

For each endpoint q of an arbitrary Li, we define 

' $kj ' Bk, if q 4- { a i b} and q is the left endpoint of an interval Jk 

i s lr'Bk, ifq 4- { a , b} and q is the right endpoint of an interval 

L ., ( 4 - 15 ) 
K°)> ifq = a 

?(1), ifq = b. 

If, for all Lj = [q[.q r ], min(|A(g/)|, |A(<3v)|) > Anpi and X(qi) ■ X(q r ) > 0, then If~\& contains 
no root of f. Otherwise, we have p < p™ ax . 



Proof If (0, 1 ) holds, then \ \g(Q) \<\g{t)\<\ \g(0) | for all t £ [0, 1] according to Lemma 7. 
It follows that 

|//(0I > |*(0I - \g(f)-fi(f)\ > ~|*(0)|- \g(t)-fl(t)\ > - («+ 1)M > 0, 

hence, / has no root in 7. Now suppose that |g(0)| < 8n/i. Since / is not contained in any J^, 
there exists a t £ [0,1] with x = a + t(b-a) £ ffl and |/(jc)| = |//(f)l < \g(t)\ + (« + < 
3^( )l + + < I6n/J,. If p > p™ ax , then from (4.13) and the definition of p™ ax , it follows 
that — log /I > p™ n = \pf + 8n — log ay] ; see the computation in (4.1 1). Hence, we have \f(x) \ < 
2~Pf . In addition, Lemma 13 (iv) and Theorem 15 guarantee that Dcm'' returns isolating intervals 
for all real roots of /, and each point in ffl has distance > a(z, , /)/(64n 3 ) from each root Zi- Thus, 
\f(x) \> (n+l )2~ p f due to Corollary 4, a contradiction. This proves (i). 

For (ii), we consider an arbitrary interval L; = [q/,q r ]. Let ?/ and t r be corresponding values in 
[0, 1] with qi = a + ti ■ w(I) and q r = a + t r - w(I). If min(|A(g/)|, |A(g r )|) > 4n/i, then 

min(|s(f z )|, \g(t r )\) > min(\X(q,)\,\X(q r )\) - (n + l)fl > 2n\i. 

Namely, for qi = a, we obviously have \g(ti) \ — \X(qi)\; otherwise, \g(ti) \ > \fi(ti)\ — (n + > 
\X(qi)\ — (n + l)jj,. For q r , an analogous argument applies. If, in addition, A (qi) -X(q r ) > 0, then 
g(t[) -g{t r ) > as well because A (qi) and X(q r ) have the same sign as g(t[) and g(t r ), respectively. 
Since we assumed that g is monotone on [0, 1], it follows that \g(t)\ > 2nfi for all t £ [f/,f r ]. This 
shows that \ fi{t) \ > \g(t)\ — (n + l)n > for all f £ [ti,t r ], thus the first part of (ii) follows. For the 
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second part, suppose that p > p™ ax . Then, B k > 2 p f > 4-njX for all k and |/(x)| > 2 p f(n + 1) 
for all x e ^ according to Corollary 4 and Theorem 15. Thus, if a <G ffl, we have 

1^(0) | > |//(0) | - (n + 1)m = \f(a) \ - (n + 1)m > 2"/ - (n + 1 )n > 4njU. 

An analogous argument applies to b. It follows that \X(q)\ > 4-njj. for all endpoints g of an ar- 
bitrary interval L; = [qi,q r ]. It remains to show that X(qi) ■ X(q r ) > 0. We have already shown 
that > 4«jU for each endpoint q, thus, f(q) must have the same sign as X(q). Namely, 

if q £ {a,b}, then f(q) differs from X(q) > 4nfi by at most (n + 1)jU < 4«jU, and, for q ^ 
{a,b}, we have sign(A(#)) = Sjy or sign(A(g)) = depending on whether q is the left or 
the right endpoint of an interval J k . Since p > p™ ax , ^ contains no root of /, thus, we must have 
X(q l )-X(q r )=f(q l )-f(q r )>0. □ 

We can now formulate the subroutine Certify'' (see Algorithm 3 in the Appendix for pseudo- 
code). Certify p is similar to DCM P in the sense that we recursively subdivide /q = (— j, j) into 
intervals / and consider corresponding p/-binary approximations // of fj. Then, in each iteration, 
we aim to apply Lemma 16 in order to certify that I<~)& contains no root of / or p < p™ ax . 
Throughout the following consideration, we assume that 

Of 

CERTIFY P never produces an interval / of width w(I) < —h^. (4.16) 

8n z 

We will prove this fact in Theorem 17 (ii). Again, we mark comments which should help to 
follow the approach by an "//" at the beginning. 



CERTIFY P . In a first step, we choose a (p + n+ l)-binary approximation / of / and evaluate 
/(— 2 +x). Then, the resulting polynomial is approximated by a (p + l)-binary approximation 

//„ € \f{-\+x)] 2 -p-u thus, //„ e [fi Q h-p according to Lemma 1. 

CERTIFY P maintains a list srf of active nodes (/,//, pi), where / = (a,b) C Iq is an interval, 
// approximates // to pj bits after the binary point and p +21ogw(/) < P/ < p. We initially start 
with srf := {(/o,// ,p)}. For each active node, we proceed as follows: 

(1) Remove (7,//,p/) from si ' . 

(2) If ln& = 0, do nothing (i.e., discard I). Otherwise, compute ^ 2 (0, 1). 

// IflC^M — 0, I is contained in one of the isolating intervals Jk, hence, we can discard I. 

(3) If t*' /2 (0,1) > -2- p '+ 2 «, check whether 

|//(0)+2- p ' +2 «| >2-P' +5 n 2 . (4.17) 
If (4.17) holds, do nothing (i.e., discard /); otherwise, return "insufficient precision". 

// For g(x) := fi(0) + 2~ pl+2 n G L//]2-P/+V tne predicate ^^ 2 (0,1) holds. From our 
assumption on w(I), we further have pi > p + 21ogw(/) > p — 2(3 + 21og« — logoy), 
and, thus, 2-P' +3 n < 2-P- 2 ( 4 "- lo g^). It follows that g fulfills the condition (4.13) from 
Lemma 16 and, therefore, I contains no root of f if (4.17) holds; otherwise, O < ap ax . 
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(4) If f^ 2 (0, 1) < -2-P' +2 n, compute h(x) = Z"= hiX l := (1 +*)"(//)' (t^) and consider the 
following distinct cases: 

(a) If hi > -n2 n ~P' for all i (or k < n2 n ~P' for all i), consider 

g( X ):=fj( X )+n2"-P'- X E [f,] n2 n- Pl 

(g(x) := fi(x) + n2 n ~ Pl ■ x, respectively). Then, for each interval L, = [qi,q r ], deter- 
mine X(qi) and X{q r ) as defined in (4.15). If min(|A(^)|, \X(q r )\) > n2 n+2 -P> and 
X(qi) -X(q r ) > for all L,, discard /; otherwise, return "insufficient precision". 

// Suppose hi > —n2 n ~ Pl for all i and g(x) := fi(x) + n2"~P'x. Then, the polynomial 
(1 +x) n {g)'{^r x ) = (1 + x) n (//)'( ) +n2 n ~P'{\ +xf has only positive coefficients. 
It follows that var(g', (0, 1)) = and, therefore g is monotone on [0,1]. In addition, 
from our assumption on w(I), we have n2 n ~P> < 2~f , ~ 2 ( 4 "~ lo S /). Hence, we can ap- 
ply Lemma 16 (ii) to g which guarantees that IC\& does not contain a root of f if 
min(|A(<7/)|, |A(g r )|) > n2 n+2 ~ pl and X(qi) -X(q r ) > for all L,- = [qi,q r ]. If one of the 
latter two inequalities does not hold, then a < (7™ ax . The case hi < n2"~ p ' for all i is 
treated in exactly the same manner. 

(b) If there exist hi and hj with hi < —n2 n ~ pl and hj > n2"~ Pl , then / is subdivided into 
// := (a,m/) and I r := (mi,b). We add {Ii,fi p pi — 1) and (I r ,fi r ,pi — 2) to si ', where 
fi t is an p/-binary approximation of //(j) and fi r an (p/ — 1) -binary approximation of 
fi{^Y~); see Step 4 (b) of Dcm p for details. If p/ < 2, return "insufficient precision". 

// Due to Lemma 1, we have fi t <G [fi^-Pi-i and fj r € \fj r ] 2 - Pl -2. Hence, by induction, it 
follows that p +21ogw(/) < p/ < p for all active nodes. 

Certify'' stops when srf becomes empty. If Certify'' returns "insufficient precision", we know 
for sure that a < a? ax . Otherwise, the region of uncertainty & contains no root of /. 



The following theorem proves that our assumption (4. 16) for the intervals produced by Certify'' 
is correct. Furthermore, we show that CERTIFY* 5 is also efficient with respect to bit complexity 
matching the worst case bound obtained for Dcm p ; see Theorem 14. 

Theorem 17. For a polynomial f as defined in (2.2) and an arbitrary p € N, 

(i) CERTIFY^ does not produce an interval I of width w(I) < Ofj (8« 2 ) and induces a recursion 
tree of size <9(£/ + ralogn). 

(ii) CERTIFY'' needs no more than 0(n{Lf + nlogn)(«r + T + p — log oy)) bit operations. 
(Hi) For a > of ax , CERTIFY P succeeds. 

Proof. An interval / is only subdivided if ty 2 (0, 1) < —2~ pl+2 n (Step (3)) or if there exist co- 
efficients k and hj ofh(x) = YH^hx 1 = (1 +x) n (//)'( ^) with k < -n2 n ~P' and hj > n2 n ~P' 
(Step 4 (b)). In the first case, we must have ty 2 (0, 1) < since |f 3 / ' 2 (0, 1) -t(' /2 (0,l)\ < 2-P' +2 n, 
hence, Ty 2 (0, 1) does not hold. For the second case, we have var((//)', (0, 1)) =^ since cor- 
responding coefficients of h(x) = (1 +*)"(//)'( jt^) and (1 +x) n (fi)'( j^) differ by at most 
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n2"-P', thus, var(/V) = var((/')/, (0, 1)) = var((//)', (0, 1)) ^ 0. Hence, the first part of (i) fol- 
lows from Lemma 9 (iv) and (iii) is then immediate from the remarks in the above description of 
Certify* 5 . Namely, Certify p only returns "insufficient precision" if p < p™ ax and guarantees 
that f(x) for all x g 2%, otherwise. For the second part of (i), we remark that, due to the above 
argument, an interval / is terminal if the disc A 2nw /[\(m(I)) does not contain a root £, of / with 
cr(£ ,/) < An 2 w{I). In [33, Section 4.2], it is shown that the recursion tree T(f') induced by the 
latter property] 4 [has size 0(Lf + nlogn). Hence, the same holds for the recursion tree induced 
by Certify'' which is a subtree of T(f'). Finally, (iii) follows in completely analogous manner 
as the result on the bit complexity for Dcm p as shown in the proof of Theorem 14. □ 

Eventually, we present our overall root isolation method RISOLATE. It applies to a polynomial 
F as given in (1.1) and returns isolating intervals for all real roots of F. 



MSOLATE: Choose a starting precision p 6 N (e.g., L = 16) and run Dcm'' on the polynomial 
/ as defined in (2.2). If Dcm p returns "insufficient precision", we double p and start over again. 
Otherwise, Dcm'' returns a list G — { (J^, Ski,Sk nBk)}k=i m with isolating intervals for some 
of the real roots of /. If Certify'' returns "insufficient precision", we double p and start over 
the entire algorithm. If Certify'' succeeds, the intervals J% = (c£,<4) isolate all real roots of /. 
Hence, we return the intervals (2 r+ c*,2 +1 rfjt), k = 1, . . . ,m, which isolate the real roots of F. 



The following theorem summarizes our results: 

Theorem 18. Let F be a polynomial as given in (1.1). Then, RISOLATE determines isolating 
intervals for all real roots of F and, for each of these intervals J containing a root E, of F, it 
holds that 

^p- <w (J)<2na^,F). 

RISOLATE demands for coefficient approximations of F to 0(Y.f +hTf) bits after the binary 
point and the total cost is bounded by 

<5(n(£ F + nY F f) = 6{n{Y. F +nl) 2 ) 

bit operations. For F G Z[jc], the bound on the bit complexity writes as (5(n 3 T 2 ). 

Proof. It remains to prove the complexity bounds and the claim on the width of the isolating 
intervals. According to Appendix 6. 1, the computation of an approximate logarithmic root bound 
r G N as defined in Section 2.2 amounts for O^nTp ) 2 ) bit operations. For a certain precision p, 
the total cost for running Dcm p and Certify p is bounded by 

6(n(E/+nlogn)(«r+T+p-logO/)) = 0(n(Lp + «r)(«r + T + p -logo*-)) 

bit operations; see Theorem 14 and Theorem 17. Since we double p in each step and succeed 
for p > pf ax , p is always bounded by 2pJ? ax = 0(£/+n) = 0{L F + nV). It follows that the 
total costs are dominated by the cost for the last run which is 0(n(Lf + nF)(nF + T + E^ )). 

4 In [33, Section 4.2], T(f') is defined as subdivision tree obtained by recursive bisection of the interval (— j, |) in 
accordance with the following rule: At depth h G No, an interval / = (— \ + r'2~''~' , — X + (/+ l)2~ h ~ l ) is subdivided if 
and only if var(/',7) ^ and A 2 s n 5 H ,( ; ) (»*(/)) contains a root | of / with separation o"(^,/) < 2 7 n 5 »'(7). For the given 
situation, we can alternatively define T(f') as the (even smaller) tree obtained by recursive bisection of (— i, i), where 
an interval I is subdivided if var(/,/') ^ and 6q MW (n{m{T)) contains a root § of / with a(^,f) < 4n 2 w(I). 
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Furthermore, we have to approximate the coefficients of / to 0(Lf +n) = 0{Lf + nF) bits after 
the binary point. Hence, the coefficients of F have to be approximated to 0(T.F+nF + T) bits after 
the binary point; see Section 2.3 for more details. From our construction of / and F, it holds that 
r < 41ogn + Tf and % = 0(nFf) (see Appendix 6.1), hence, we can replace F by Ff and further 
omit T in the above complexity bounds. For the special case where F is an integer polynomial, 
the bound on the bit complexity follows from = 0{m); see Appendix 6.2. The estimate on the 
size of the isolating intervals is due to the following consideration: An interval / which contains 
the root z = ^frr of / is not subdivided by Dcm'' if w(I) < a(z,/)/(8n 2 ). Hence, any interval 
which is returned by Dcm'' as an isolating interval for z is the extension I = (a — ^5^,^ + 
of an interval (a,b) with w(I) > a(z,/)/(16n 2 ), thus w(J) = 2 T+x w(I) > o(%,F)/(16n 2 ). 
From our construction, the ^^-neighborhood of / isolates z as well and, thus, w(J) = 2 r+1 w(I) < 
4Fno{zi,f)=2no(%,F). □ 

4.4. Some Remarks 

4.4.1. On the Complexity Analysis for Integer Polynomials 

We remark that in order to achieve the complexity bound (5(n 3 T 2 ) for integer polynomials, 
the subroutine CERTIFY^ and its analysis is not needed. Namely, due to our considerations in 
Appendix 6.2, we can compute upper bounds for Lf (in terms of n and t) and, thus, also an 
upper bound p*(n, t) for p™ ax which matches p™ ax at least with respect to worst case complexity. 

Then, according to Theorem 15, it is guaranteed that Dcm p *(" ^ computes isolating intervals for 
all real roots of /. Unfortunately, this approach cannot be considered practical at all because such 
upper bounds usually tend to be much larger than the actual p™ ax . We would like to emphasize 
on the fact that our algorithm is output sensitive in the way that it demands for a precision which 
is not much larger than Pf, hence, our algorithm chooses an almost optimal precision. Without 
giving an exact mathematical proof, we conjecture that, for any bisection method, the bound on 
the bit complexity as achieved by our algorithm is optimal (up to log-factors). Namely, the bound 
<9(«t) on the precision as well as the bound 0{nx) on the size of the recursion tree cannot be 
lowered for Mignotte polynomials. 

4.4.2. On Efficient Implementation 

We formulated our algorithm in a way to make it accessible to the complexity analysis but still 
feasible and efficient for an implementation. Nevertheless, we recommend to consider a slight 
modification of our algorithm when actually implementing it. 

For our certification step CERTIFY, the most obvious modification is to only subdivide the 
region 8% instead of the entire interval (— j, 5). More precisely, 3% decomposes into intervals Lj 
"in between" the isolating intervals J^. Then, we approximate the polynomials fa. to p bits after 
the binary point and recursively proceed each Lj in a similar way as proposed in Certify p . 
An experimental implementation of our algorithm in MAPLE has shown that following this ap- 
proach the running time for the certification step is almost negligible whereas, for the original 
formulation, it is approximately of the same magnitude as the running time for DCM P . 

Furthermore, we propose to also use the inclusion predicate based on Descartes' Rule of 
Signs. With respect to complexity, our inclusion predicate based on the 3^', 2 -test (see Corol- 
lary 8) is comparable to Descartes' Rule of Signs, where we check whether / has exactly one 
sign variation for a certain interval. However, in practice, this subtle difference is crucial because 
already log n bisection steps more for each root may render an algorithm inefficient. As an alter- 
native, for each interval / in the subdivision process, we propose to check whether there exists a 
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suitable p/-approximation g of // with var(g, (0, 1)) = 1. Namely, if there exists such a g, then we 
can proceed with fi : = g which has exactly one root in /. Thus, it is easy to refine / (via simple 
bisection or quadratic interval refinement) such that ,5^ 2 (0,2) holds as well. 

We finally report on an interesting behavior of the proposed method. It is easy to see that, 
for small intervals /= (a,b), the leading coefficients of fi(x) = f(a + w(I)x) are considerably 
smaller than the first-order coefficients. Since we only consider a certain number p/ < p of bits 
after the binary point, the approximations // are usually of lower degree than fj. As a conse- 
quence, the cost at such an interval is tremendously reduced because we have to compute the 
polynomial //+ = //(^ + (1 + which is expensive for large degrees. In particular, for a 
polynomial with two very nearby roots (such as Mignotte polynomials), this behavior can be 
clearly observed. More precisely, when refining an interval / which contains two nearby roots, 
the degree of // decreases in each bisection step and eventually equals 2 for / small enough. We 
consider this behavior as quite natural because // implicitly captures the information on the loca- 
tion of the roots in a neighborhood of / whereas the influence of all other roots becomes almost 
negligible. In comparison to previous methods, the proposed algorithm exploits this fact. 

5. Conclusion 

We presented a new complete and deterministic algorithm to isolate the real roots of an ar- 
bitrary square-free polynomial F with real coefficients. Our analysis shows that the hardness of 
isolating the real roots exclusively depends on the location of the roots and not on the coefficient 
type. Furthermore, the overall running time is significantly reduced by considering approxima- 
tions at each node of the recursion tree. In particular, for integer polynomials, we achieve an 
improvement with respect to worst case bit complexity by a factor n = degF compared to the 
best bounds known for other practical methods such as the Descartes or the continued fraction 
method. The latter is due to the fact that exact arithmetic produces too much information for the 
task of root isolation and, thus, a significant overhead of computation. We remark that recent 
work [20] shows corresponding results for the task of further refining given isolating intervals. 
Since we are aiming for a practical method, we formulated our algorithm in the spirit of the 
Descartes method. We are convinced that because of its similarities to the latter exact method 
and because of its usage of approximate and less expensive computations, it will prove to be 
efficient in practice as well. We plan to implement our algorithm to verify this claim. A first 
promising step [37] has already been made. Finally, univariate root isolation constitutes an im- 
portant substep in cad (cylindrical algebraic decomposition) computations. In combination with 
a recent result on the complexity for real root approximation [20], it is possible to obtain a bound 
on the worst case bit complexity of topology computation of planar algebraic curves [21] which 
crucially improves upon the existing record bounds from [9,19]. 
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6. Appendix 



6.1. Approximating Tp 

In this section, we show how to compute an integer approximation T e N of the exact loga- 
rithmic root bound Tp := log(max; |^, |) with Tp < T < 41ogn + Tp. We further prove that this 
computation can be done with 0((nTp) 2 ) bit operations. As a byproduct, x = 0(nTp). 

Consider the Cauchy polynomial F c (x) := \A n \x n — £"=o \Ai\x 1 of F. Then, F c has a unique 
positive real root and it holds that 

2 r> <( e* < _ n _. 2 r F <2n-2 r "; 
ln2 

see [12, Proposition 2.51]. Thus, it follows that F c (x) > for all x > ii2 Yf+1 and, in particular, 
\A n \{n2 TF+l f > \A i \{n2 TF+l ) i for all i. From the definition of X on page 2, we have x = 1 or 
there exists an ;'o with \A io \/\A n \ > 2 T_1 . The first case is trivial and, in the second case, we have 
(2 r >+ 1 n)"- ! o > 2 T_1 . Thus, x = 0(nT». 

We now aim to compute an approximation of E, * via evaluating F c at x = 2,4, 8, Then, for 

the smallest k £ N (denoted by ko) with F c (2 k ) > 0, we must have I;* < 2 k ° <2^*. However, 
since F has approximate coefficients, these evaluations cannot be done exactly. The idea is to use 
interval arithmetic with a certain precision p (fixed point arithmetic) such that 25 (F c (2 k ),p) < 1 , 
where *B(E,p) is the interval obtained by evaluating of a polynomial expression E via interval 
arithmetic with precision 2~ p for the basic arithmetic operations; see [20, Section 4] for details. 
We initially start with k = 1. If <B(F c (2 k ),p) contains zero or A < for all X e ( B(F c (2 k ),p), 
we proceed with k+ 1. Otherwise, we must have feo < A: < fco + 1: The left inequality is obvious. 
For the right inequality, we remark that 2 k ° +l has distance larger than 1 to all roots of F c and, 
thus, F c {2 k « +1 ) > 1. Hence, X > for all X e < B(F c (2 A: o +1 ),p). It follows that T:=k fulfills 

T F < r < k + 1< log(4<in < log(8«2 r ' r ) < 41ogn + r F . 

It remains to bound the cost for the interval computations of < B(F C (2 k ),p). Since F c has 
coefficients of size less than 2 T , we have to choose p such that 

2-P +2 (n+l) 2 2 T+nk < 1 

in order to guarantee that w( < B(F c (2 k ),p)) < 1; see [20, Lemma 3]. Then, p is bounded by 0{x + 
nk) and, thus, each interval evaluation needs 0(n(x + nk)) bit operations. From our above con- 
siderations, we have k < ko + 1 = 0(\ogn + Ff), hence, the total cost is bounded by 0{nTf{x + 
nTp)) = d((nrf) 2 ) and we need approximations of F to 0(x + nFf) = O'liTp) bits after the 
binary point. 

6.2. Integer Polynomials 

For an integer polynomial F G Z[x] as given in (1.1), we aim to show that Lp = 0{nx). We 
proceed in two steps: First, we cluster the roots of F into subsets consisting of nearby roots. 
Second, we apply the generalized Davenport-Mahler bound [10,12] to the roots of F. Eventually, 
the above result follows. 

W.l.o.g., we can assume that o(£i,F) <...,< o(% n ,F). For h e N, we denote i(h) the maxi- 
mal index ;' with a(§ ( -,F) < 2~ h and R = R(h) := } the corresponding set of roots ^ 
with a{E,i,F) < 2~ h . lfh< log(l /dp), then R contains at least two roots. We are interested in a 
partition of R into disjoint subsets R\,...,Ri that consist of nearby points, only. 
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Lemma 19. Suppose that h < log(l/c>). Then, there exists a partition ofR := R(h) into disjoint 
sets Ri,...,R[ such that \Rj\ > 2 for all i e {1,...,/} and \£, < n2~ h for all %, %' e R ( . 

Proof We initially set R\ := Then, we add all roots § to #1 that satisfy |§ - < 2~ /l . 
For each root in R\, we proceed in the same way. More precisely, for each £ e R\, we add those 
roots 4' e to /?i with || — <i;'| < 2 _A . If no further root can be added to R\, we consider the 
set R\R\ of the remaining roots and treat it in exactly the same manner. Finally, we end up with 
a partition R\,...,Ri of R such that, for any two points in any their distance is less than or 
equal to (\Rj\ — l)2~ h < n2~ h . Furthermore, each of the sets Rj must contain at least two roots as 
<j(£i,F) <2~ h for all i = l,...,i(h). □ 

We now consider a directed graph 5^- on each R{ which connects consecutive roots of Ri in 
ascending order of their absolute values. We define <S := (R,E) as the union of all S^-. Then, <S is 
a directed graph on R with the following properties: 

(1) each edge (a,j3) e E satisfies \a\ < |J3|, 

(2) Sf is acyclic, and 

(3) the in-degree of any node is at most 1 . 

Hence, we can apply the generalized Davenport-Mahler bound [10,12] to ( S: 

1 { vA \\ n ' 2 



As each set R{ contains at least 2 roots, we must have i(h) > #E > i(h) /2. Furthermore, for each 
edge (a,P) € E, we have \a - J3| < «2~\ It follows that 

i(A) 



( " j 2 > (V«+T2 T )«- 1 ' ^ « J "W > (" + 1)"2" T 'U 2 / 
and, thus, 

2n(T + log(n+l)) 2n(T + log(n + l)) 

i(h) < < . 

log3+log« + /z /z 

It directly follows that log ( 1 jap ) < n ( T + log (n + 1 ) ) + 1 since, otherwise, there would exist an 

h with n(r + log(n + 1)) < h < log(l/o>) and i(h) < 2 which is not possible. For the bound on 

, it suffices to consider only the roots ^1, . . . , ^ with separation < 1/2 since all other roots 

contribute with at most n to the sum Lp. Since 

k [«(T+log(«+l))] [n(T+log(n+l))l j 

-J>g<r($,F)< £ /(/z)<2«(T + log(n + l)) £ - = 0(«Tlog(«T)), 

i=l A=l /i=l 

it follows that Z f = 0{nx). 
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6.3. Algorithms 



Algorithm 1 Dcm 

Require: polynomial / = Y*o<i<n a i x ' € as defined in (2.2) 

Ensure: returns a list & of disjoint isolating intervals for all real roots of / 

{only in the REAL-RAM model} 

h-={-\,\) 
f k {x):=f{-\+x) 

srf := { (Io,fi ) }; ^:=0 {list of active and isolating intervals} 

repeat 

(/,//) some element in srf with / = (a,b); delete (/,//) from srf 

//+= = //(-£ + (l + Si) *) and fh W = E?M) V == (1 +*)" ■ //+ (l^) 

if var(// + ) =0then 

do nothing 
else 

if f^' (0,2) >0then 

*:=sign/ /+ (0)-/ /+ (l) 
if s > then 
do nothing 
else 

if I + does not intersect any interval in & then 

add /+ to G 
else 

do nothing 
end if 
end if 
else 

subdivide / into // := (a, mi) and I r := (mj,b) 
fi, :=// (f) and fa :=/, (^) = / /( (x + 1) 
add (/,,//,) and (I r J, r ) to ^ 
end if 
end if 
until is empty 
return ^ 
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Algorithm 2 Dcm p 

Require: polynomial / = Y,o<i<n a i x ' € M.[x] as in (2.2) and a p e N 

Ensure: returns "insufficient precision" or a list & = {Jk,Sk,i, s k,nBk} of disjoint isolating in- 
tervals Jk = {c^.dk) for some of the real roots of / (and Sk ; = signf(cfc), s^ r = sign/(<4) and 
0<B fe <min(|/(c ft )|, 



h'-={-\,\) 

fa (p + n + 1) -binary approximation of / 

fi a (p + l)-binary approximation of /(— j +x) 

^:={(k,f k ,p)};G:=& 



// G [// ] 2 -p } 
{list of active and isolating intervals} 



repeat 

(/,//, pi), where /:= (a,b), some element in srf ; delete (/,//, pi) from 
/ /+ (x) :=/, (-£ + (1 + and fc» = E» =0 %<' := (1 +*)»■ j* + (^) 
if hi > -2 n+2 -Pi for all i or A,- < 2 n+2 -P< for all i then 

do nothing 
else . 

if > -n2 n+l -P< then 

fi{x):=fi{x)+n2 n+l -P< -x 

X- :=f I+ (0)-2"- l -P>, A+ :=/ /+ (l) + (4«+ 1)2"- 1 -'" and A :=f,(-l/n) -2 n+1 ~P>. 

if /= (a — m^-,b+ intersects no interval 7 for all {J ,sj j,sj^ r ,Bj) S ^ and A~ • 
A+ < and rnin(|A-|, |A+|) > n2 n+3 -P> and |A| > „2 2 deg(/ / )+7+„- P; then 
add (/, sign(A - ) , sign(A+) , rnin(| A" | , | A+ 1 ) - 2 n+3 -P'n) to G 

{=>■ I contains a root t, of / and the -neighborhood of / is isolating for <^ } 

else 

do nothing {/is already isolating for E, } 

end if 
else 

do nothing 
end if 
else 

if pi < then 

return "insufficient precision" 
else 

if pi < 2 then 

return "insufficient precision" 
else 

Subdivide I into 1/ := (a, mi) and I r := (mi,b) 

fii an p/-binary approximation of // ( § ) {^^tf/fU-i)} 
// r an (p/ - 1 ) -binary approximation of // ( ^ ) { => f h e [//,. ] 2 _ (P; _ 2) } 

Add (/, J h ,pi-\) and (/„ fi r , pi -2) to 
end if 
end if 
end if 
until srf is empty 
return 6 
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Algorithm 3 Certify* 5 

Require: polynomial / = Y.o<i<n a i x ' <= as defined in (2.2), an p € N and the list 6 = 

{(Jk,Sk,l,Sk,r,B k )} k =i,...,s returned by Dcm p . 

Ensure: returns "insufficient precision" or the list Jzf — {Jk}k=i,...,s of isolating intervals with 
the guarantee that, for each real root of /, there exists a corresponding interval in Jz? . 

f an (p +n + 1) -binary approximation of / 

fh an (P + 1) -binary approximation of /(- \+x) f h) e [fi ] 2 -p } 

£/:—{ (Io,fi ,p) } {list of active intervals} 

repeat 

(/,//, Pi), where /:= (a,b), some element in a/; delete (/,//, pi) from si. 

if/n^' = Ui=i z 'i = 0then 

do nothing 
else 

if ^ 2 (0,1) > -(n + l)2-P> +1 then 

if |//(0) +2-P' +2 n\ > (n+l)2-P'+ 5 then 

do nothing {7 contains no root of /} 

else 

return "insufficient precision" {p < p™^} 

end if 
else 

h(x):=U =0 h i x i = (l+xr(fi)'(^ x ) 

if hi < n2"-P> for all i (or h t > -n2 n ~P' for all i) then 

g(x) := fi(x) - n2"-P> (or g(x) := fi(x) + n2 n ~P', respectively); 

if foreachL, = [qi,q r },mm(\X(qi)\,\X(q r )\) >n2 n+2 ~P' taaA X{qi) ■ k{q r ) <0then 
do nothing {7n,^ contains no root of /; X(qi), X(q r ) defined as in (4.15)} 

else 

return "insufficient precision" {p < pjj™} 

end if 
else 

if pi < 2 then 

return "insufficient precision" 
else 

Subdivide / into // := (a, mi) and I r := (mj,b) 

fi, an p/-binary approximation of // (§) //, e [// ; ] 2 -(p/-i) } 

/>, an (p/ - l)-binary approximation of // (^) f h e Lfr] 2 -(p/-2)} 

Add {Ii,fi p pi- 1) and (I r ,fl,pi-2) to 
end if 
end if 
end if 
end if 
until is empty 

return "certification successful" {The region of uncertainty ^ contains no root of /} 
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